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This paper investigates the dynamical generation of entanglement in scattering systems, in par- 
ticular two spin systems that interact via rotationally-invariant scattering. The spin degrees of 
freedom of the in-states are assumed to be in unentangled, pure states, as defined by the entropy of 
entanglement. Because of the restriction of rotationally-symmetric interactions, perfectly-entangling 
S-matrices, i.e. those that lead to a maximally entangled out-state, only exist for a certain class of 
separable in-states. Using Clebsch-Gordan coefficients for the rotation group, the scattering phases 
that determine the S-matrix are determined for the case of spin systems with a = 1/2, 1, and 3/2. 
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I. INTRODUCTION 



(f-) ■ The study of entanglement is of central importance to quantum information theory. Entanglement is the resource 
\ for many of its proposed applications, such as quantum computation and quantum teleportation A variety 
C*~) i of experimental procedures for generating entangled quantum systems have been demonstrated and/or proposed and 
the search continues for new methods that will be stable, scalable, and efficient. To guide and explain this work, it 
is important to have theoretical studies of how entanglement is generated in quantum states by the application of 
dynamical operators. 

, This work investigates the dynamical generation of entanglement in scattering spin systems by considering the 
S-matrix operator and how global symmetry properties restrict its form. This is an complement and alternative to 
two related methods for investigating the generation of entanglement by dynamics. One can consider the general 
properties of global unitary operators by looking at classes of unitary operators in the state space of the quantum 
system of interest, typically the tensor product of qubits. See |3j for work related to the approach of this article 
and for a comprehensive review and bibliography. In another approach, one can investigate possible interaction 
Hamiltonians and their exponentiations. These Hamiltonians could be purely internal or could incorporate the action 
of outside fields. The papers Q and provide good examples of this approach, and many more examples exist which 
O^ 1 are specific to particular experimental configuration. 

The approach of this article specifically explores the generation of entanglement in the non-relativistic, elastic scat- 
tering of two distinguishable particles by a central force interaction. Non-relativistic two-body interactions dominate 
the dynamics of a gas of trapped ultra-cold atoms, for example, and applications of quantum information theory to that 
system 0, H> require an understanding of dynamical entanglement by scattering 0] • More generally, these results 
apply to any two particle "scattering-like" experiment, i.e., a bi-partite spin system where a spherically-symmetric 
interaction between the two spins can be turned on and off (see example in |TlJ). Systems that are asymptotically non- 
interacting can be cast into the form of a scattering problem and treated with the techniques below. Such a sequence 
could be arranged via controlled interactions, but also appears naturally in the case of finite-range interactions. 

In general, the description of entanglement in scattering systems requires that one considers entangled states of 
with continuous degrees of freedom (see |12| and references therein, and the review |l3)p. However, for non-relativistic 
particles with central interactions, there is no mixing between orbital and intrinsic angular momentum [l4j |. Then the 
S-matrix can be decomposed into partial waves S s i(p, E) labeled by total spin s and by orbital angular momentum I 
and the S-matrix is diagonal in total momentum p and energy E. The scattering dynamics do not mix partial waves 
with different values of s and I, and so within each partial wave of orbital angular momentum the entanglement of 
the spin degrees of freedom is separable from the translational degrees ^5|. As a result, this article only considers 
entanglement in the spin degrees of freedom. One cannot make this separation between the translational and rotational 
degrees of freedom if the interactions are non-central 0] or if the system is relativistic 0, . 
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In the case of bi-partite quantum states with finite degrees of freedom, the quantification of the entanglement in 
a pure state is unambiguous. Then one may ask, whether there exist S-matrices that are perfect entanglers, i.e. an 
operator that takes an initially separable state into a maximally entangled state Q- The main result proved here is 
that rotationally- invariant, perfectly-entangling S-matrices only exist for a particular type of separable initial state. 
When they do exist, the scattering phase shifts that determine the S-matrix in the spin sector can be explicitly 
calculated. For scattering spin systems like those described above, this paper calculates the in-states and phases 
necessary for maximal entanglement of spin systems with spins of a — 1/2, 1, and 3/2. Reversing the idea, this shows 
how information about scattering interactions can be gained by looking at the entanglement of the out-state within 
particular partial waves of orbital angular momentum. 



II. DYNAMICAL ENTANGLEMENT 



Consider two quantum systems with the same finite- number of levels d with Hilbert space Ti.^ — Ti-d <8> T~td- The 
entropy of entanglement for a pure state G TL d 2 is 

Ety) = S{pi) = S( P2 ) (1) 

where p± = tr2[| ?/')('</' I] is the density matrix for system 1 that remains after a partial trace over system 2, and 
S(p) = — tr[plogp] is the Von Neumann entropy of the density matrix p. Conventionally, the logarithm in the entropy 
is taken in base 2, but for our purposes it is better if it is taken base d. Then the entanglement is bounded by 
< E(ijj) < 1. A pure state of the form 

\4>) = \</>i)®\<h) 

= E a ti) x E « ( 2 ) 

3=0 fc=0 



is unentangled. The reduced density matrix is pi — \4>i){4>i\ and E(<j)) = 0. States of the system with the form 

1 

Vd 



I (d-l) 

W) = —, E e<aj li>®K>> ( 3 ) 

3=0 



where nj is the j-th element of permutation 7r 6 Sd of the numbers {1, d} and <Xj G R, have maximum entanglement 
E(ip) = logoff = 1. These states have reduced density matrices p\ = p 2 — (l/d)I<j. 

How can a maximally-entangled state © evolve from a unentangled pure state J5J? In other words, what operators 
are perfect entanglers for such a system. Mathematically, any unit-normalized state in H d 2 can be transformed to 
any other unit-normalized state in H. d 2 by unitary transformation U G \J(d 2 ). Since a global phase in not physically 
meaningful, it is sufficient to consider SU(d 2 ). The subset of local operators U(c?) x U(d) cannot change the entropy 
of entanglement of a pure state, so perfect entanglers must be elements of the set P d = SU(d 2 )/U(rf) x XJ(d). Zhang 
et al. consider the system of two qubits and find that half of P 2 are perfect entanglers. 

From another perspective, since U(iV) is a connected matrix Lie group, it is possible to express every U G U(d 2 ) as 



U = exp(iHiti) exp(iEf 2 *2)"- exp(iH m t m ) (4) 

for some finite number of d 2 x d 2 Hermitian matrices {Hi, H2, ...H m } [Is| . In principle, one could imagine some series 
of interaction Hamiltonians, switched on and off at certain times, that when exponentiated lead to time evolution 
operators that are perfect entanglers [3j. However, in a given physical system, symmetries may limit what kind of 
Hamiltonians can appear in or equivalently many U G U(d 2 ). 

The example we consider in depth here is scattering by central forces of two systems with spin a. Assume that 
in the limit t — > ±00, the two systems are not interacting. Then one can define the unitary scattering operator, the 
S-matrix, that transforms the unentangled in-state \<f>) to the maximally-entangled, out-state \<fi'}: 

m=s\ct>). (5) 

Scattering interactions are spherically symmetric so the S-matrix must commute with the total angular momentum 
angular momentum operator. Since it is a central interaction, it must also separately commute with the total spin and 
orbital angular momentum operators. The spin sector of the interaction can be separated from the other degrees of 
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freedom and will therefore be identified with Tid 2 , where d = 2a + 1, and the rotations R € SO (3) will be represented 
by the tensor product of the single-particle spin representations D a (R) x D a (R). As a result, the set of possible 
perfectly-entangling S-matrices is much reduced by the symmetry: they must global, unitary matrices that commutes 
with every D a (R) x D a (R). Additionally, it will be shown that only for a certain set of in-states \<fi) do S-matrices 
exist that are perfect entanglers and the possibilities can be enumerated. 



III. USEFUL BASES FOR ANALYZING SPIN SYSTEM DYNAMICAL ENTANGLEMENT 

There are two bases that will be used for the states in 7^2 . The first is the direct product basis denoted by either 
|Mi>M2) or |xi)X2)- These are the eigenvectors of individual angular momentum 3-component operators: 

s 3 Vl>M2) = Ml|A*l,M2) 

S 3 2 Vl.M2) = M2|Ml,M2)- (6) 

The direct product basis is useful because it is separable like the initially-unentangled state and because measurements 
are local operations. Also, the entropy of entanglement requires the partial trace, which is straightforward to evaluate 
in this basis. 

The second useful basis is the direct sum basis |sm), which are the eigenvectors of the total angular momentum 
component operators E3 and total angular momentum squared operator S 2 : 

Y>3\sm) — m\sm) 

S 2 |sm) = s(s + 1) \sm). (7) 

This is called direct sum basis because it arises in the Clebsch-Gordon decomposition of products of irreducible 
representations (IRs) of the SU(2) into a direct sum of IRs It is useful because if the interaction is spherically- 
symmetric (hence SU(2)-invariant), then the S-matrix has the form 

(s'm'\S\sm) = e 2lSa 5 ss >5 mm > . (8) 

The convention of calling the s-dependent phase 2S S comes from scattering theory and the form (jSJ can be seen as the 
consequence of the Wigner-Eckhart theorem for scalar operators or more generally as the consequence of Schur's 
lemma |20) . In the case of scattering non-relativistic particles,<L also depends on the total angular momentum, orbital 
angular momentum, and magnitude of relative momentum [15| . The direct product basis and direct sum basis are 
connected by the Clebsch-Gordan coefficients (CGCs) for SU(2), (sm|^ 1 ,/i 2 ), which can be chosen as real. 

So then the question becomes: for what states <fi of the form J5J will S<f> be a maximally entangled state of the form 
(|3J) and what must the phases S s be for this to occur? To answer this question, the state <j)' — S<j) is expressed in the 
direct product basis: 

W) = s X! a ;<A'2 iMiiAta) 

~ ^ C Xl,Xl\Xl,X2) (9) 



/ j Xl.X2 

X1X2 



where 



c xiX2 = a Mi^2^(xi,X2;Mi,M2) (10) 
P1M2 

and 

KXiiXaJMijMa) = ^2(s m \^i, (xi, X2\s m)e 2l5s . (11) 

s in 

Using this notation, the reduced density matrix for particle 1, pi = tr2(|0')(0'|), is 

Pi= EE c xx,x 2 c xl)X2 |xi><Xi|. (12) 

Xl Xi X2 

Since p\ = (l/e?)!^ maximizes S(p\), the state <f>' = S<f> will be a maximally entangled state if 

E c xi.X2 c x ' 1 , X2 = ^xi.x'r ( 13 ) 
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IV. FINDING THE IN-STATES AND PHASE SHIFTS 



The equation H13|) must be solved to find the possible phase shifts and in-states that lead to maximally- entangled 
out-states. The coefficients c(xi,X2) m depend on three things: the initial state <p through a M1 and bu 2 , the 
CGCs for SU(2), and the scattering phases S s . 

The only states for which a perfectly-entangling S-matrix exists, as will be shown below, are states of the form 

\<t>(u,\)) = U(u)\\,-\), (14) 

where U(u) = U\(u) <g> U%[u) is the direct product representation of the rotation group with u e SU(2) and A e 
{cr, a — 1, ... — cr}. Such states are the zero eigenvectors of the total spin operator E3 = (R(u)H) 3 , where R(u) is the 
image of u under the standard homomorphism SU(2) — > SO(3). 
Because [S, U(u)] — for all u £ SU(2), we have 



= U(u)S\X,-X) 



e 2i5 °(s 



m = 0|A,-A)|sm = 0) 



[/(n)^c Xl , X2 (A)| Xl ,X2), (15) 



X1X2 



where we define 

ZxuxiW = 5]e 2 ^( S m-0|A,-A)(xi,X2|sm-0) 



Then l|12ll becomes 



Pi = U 1 (u) ^| ffxi (A)| 2 | X i)(xi|^ 1 t H. 



(17) 



This is a diagonal matrix and |<7 Xl (A)| 2 are the Schmidt coefficients for the state after the interaction. If all the 8 S 
are the same phase, then |<? Xl (A)| 2 = ^xi,a because the basis transformation given by Clebsch-Gordan decomposition 
is unitary. In this case, the dynamics just evolve the in-state by a total phase and there is no entanglement. The 
reduced density matrix p\ in (|17|) will be of the maximally entangled form for a given A if and only if 

| 5xi (A)| 2 = l/dforall Xl . (18) 

If (|18|) is satisfied, then the density matrix will be a scalar multiple of the identity and commute with all rotations 
u E SU(2). It can be shown that (|17fl would not be diagonal for any other initial condition besides one of the form 
4>(u, A). Only eigenvectors of the total angular momentum component (in any direction) with m = lead to a reduced 
density matrix of the form l|17f) because only those states can have non-zero Clebsch-Gordan coefficients with every s 
from zero to 2a. Also, this shows that all maximally-entangled states that emerge from a scattering experiment will 
have the form in (|15f) . 

When the in-state kt(j)(u,X) has the form iQ , then the set of <5 S which satisfy of |g Xl (A)| 2 = 1/d for all xi <= 
{cr, cr — 1, ... — cr}, will determine a perfectly-entangling S-matrix. Explicit solutions have been found for a = (trivial), 
1/2, 1, and 3/2 and it has been found that the solutions are in fact independent of A and only depend on a. The results 
are summarized below. The phases are set so Sq = and all other phases are relative to this and all i5 s S (— tt, 7t]. 

• For cr = 0, spin entanglement is not meaningful. 

• For cr = 1/2, l|18l) for all A and xi leads to two independent equations: 



= cos 2 S\ 
= sin 2 Si . 



By fixing the global phase so that <5o = 0, the geometric structure of possible S-matrices [Si] is isomorphic to 
the circle S 1 and the four perfectly-entangling S-matrices are the points Si = ±37r/4 or Si — ±7r/4. 
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• For cr = l, (|18f) for all A and xi leads to three independent equations: 

I = ^-{7 + 6cos(26 1 )+3cos(26 1 - 26 2 )+cos{26 2 )) 
o 18 

3 = 9 Sm h 

\ = -^-(7-6cos(26i)-3cos(2(5i-2(52) + 2cos(2(52)). 
3 18 

With fixed So, the geometric structure of possible S-matrices [5i,<5 2 ] is isomorphic to the torus T 2 = S 1 x S 1 
and the eight perfectly-entangling S-matrices are the points {61,62} = {tt/12 ± ir /4, —ir/6 ± ir/2}, {—ir/12 ± 
tt/4, 7r/6±7r/2}. 

• For a = 3/2, l|18|l leads to four independent equations: 

I = (132 + 90 cos(2<5i) + 90 cos(2(5i — 26 2 ) + 18 cos(2<$i - 26 3 ) 
+50cos(2(5 2 ) + 10cos(2<S 2 - 26 3 ) + 10cos(2£ 3 )) 

\ = 450 ( 132 ~ 90cos ( 2,5 i) -90 008(2(5! - 2<5 2 ) + 18cos(2<Si - 26 3 ) 

+50cos(2(5 2 ) - 10cos(2<5 2 - 2<5 3 ) - 10cos(2£ 3 )) 
I = ^q(68 - 30cos(2 ( Ji) + 30cos(2(5i-2(5 2 )-18cos(2 ( 5i-253) 

-50cos(2(5 2 ) - 30cos(2<S 2 - 26 3 ) + 30cos(2£ 3 )) 
\ = ^( 68 + 30cos ( 2(5 i)" 30cos ( 2,5 i" 2 ' 5 2)-18cos(2 ( 5i-253) 

-50cos(2(5 2 ) + 30cos(2<5 2 - 26 3 ) - 30cos(2^ 3 )) 

With fixed 60, the geometric structure of possible S-matrices [61, 62, 63] is isomorphic to the 3-torus T 3 = 
S 1 x S 1 x S 1 and there are four sets of solutions for any value of $i E (— w, tt]: 6 2 = ±tt/2 and 6 3 — 6\ or 
62 = ±7r/2 and 6 3 = 5i ± 7r (whichever one of <5i ± 7r G (— tt, tt]). 

Solutions for a > 2 have not yet been investigated. 



V. CONCLUSION 



In summary, it has been shown that a rotational-symmetric interaction acts as a constraint on dynamical generation 
of entanglement in scattering systems. Only for certain initial states do there exist perfectly-entangling S-matrices, 
only certain phase shifts allow for such a transformation, and only a subset of maximally entangled vectors emerge as 
out-states. The implications for production of entangled states by scattering-type interactions are worth exploration 
and could guide experimentalists in constructing or searching for suitable systems and for tuning the interactions 
of those systems where possible. Additionally, this idea can be reversed, and as in classic partial wave analysis, 
en tang lement correlations could be used to find information on the phase shifts. This idea has been partially explored 
in |l2( for translational entanglement in scattering, but much work remains to be done. 

Finally, this paper can also be thought of as showing how interaction symmetries limit the possible unitary trans- 
formations and therefore limit the maximal entanglement possible for a given initial condition. Particles are elements 
of the representation spaces of space-time symmetry groups. When direct products of these representation spaces are 
decomposed into a direct sum of irreducible representations, one can represent the S-matrix in a basis that makes 
global symmetries explicit. In this case, it has been shown that for unitary irreducible representations of the Galilean 
group, the spin degrees of freedom and be separated from the translational degrees of freedom and the consequences 
of global rotational symmetry can be explored. Extensions to other space-time and interaction symmetry groups will 
be considered in the future. 



6 



Acknowledgments 

The author would like to thank Mark Byrd for fruitful discussions. 



[1] P. Shor, in Proc. 35th Ann. Symp. Found, of Comp. Sci., (IEEE, Los Alamitos, CA, 1994), p. 124. 

[2] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 (1993). 

[3] J. Zhang, J. Vala, S. Sastry, K.B. Whaley, Phys. Rev. A 67, 042313 (2003). 

[4] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, (Cambridge U. Press, Cambridge, 
2000). 

[5] W. Diir, G. Vidal, J. I. Cirac, N. Linden, S. Popescu, Phys. Rev. Lett. 87, 137901 (2001). 
[6] K. Hammerer, G. Vidal, J.I. Cirac, Phys. Rev. A 66, 062321 (2002). 

[7] D. Jaksch, H.J. Briegel, J. I. Cirac, C. W. Gardiner, and P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). 

[8] G. K. Brennen, I. H. Deutsch, and P. S. Jessen, Phys. Rev. A 61, 062309 (2000); G. K. Brennen, I. H. Deutsch, and C. J. 

Williams, Phys. Rev. A 65, 022313 (2002). 
[9] A. M. Dudarev, R. B. Diener, B. Wu, M. G. Raizen, Q Niu, Ph ys. Rev. Lett. 91, 010402 (2003). 
[10] J. Wang, C. K. Law, and M.-C. Chu, arXiv: quant-ph/0503020 

[11] C. Cohen- Tannoudji, B. Diu, and F. Laloe, Quantum Mechanics, Vol. 2, (Wiley, New York, 1977), Chap. 10, Comp. F. 
[12] A. Tal and G. Kurizki, Phys. Rev. Lett. 94, 160503 (2005). 
[13] J. Eisert and M. Plenio, Int. J. Quant. Inf., 1 479 (2003). 

[14] M. L. Goldberger and K. M. Watson, Collision Theory, (John Wiley and Sons, New York, 1964), ch. 6-7. 

[15] N. L. Harshman, Int. J. Mod. Phys. A, 20, 6220 (2005). 

[16] M. Czachor, Phys. Rev. A 55, 72 (1997). 

[17] N. L. Harshman, Phys. Rev. A, 71, 022312 (2005). 

[18] B. C. Hall, Lie Groups, Lie Algebras, and Representations, (Springer, New York, 2003), Corollary 2.31. 
[19] M. E. Rose, Elementary Theory of Angular Momentum, (Wiley, New York, 1957). 

[20] J. Fuchs and C. Schweigert, Symmetries, Lie Algebras and Representations, (Cambridge, 1997), ch. 5.4. 



